Larkin-Ovchinnikov-Fulde-Ferrell Phase in (TMTSF) 2 C10 4 Superconductor: Theory 

versus Experiment 
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We consider a formation of the Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) phase in a quasi-one- 
dimensional (Q1D) conductor in a magnetic field, parallel to its conducting chains, where we take 
into account both the paramagnetic spin-splitting and orbital destructive effects against supercon- 
ductivity. We show that, due to a relative weakness of the orbital effects in a Q ID case, the LOFF 
phase appears in (TMTSF^CICU superconductor for real values of its QfD band parameters. We 
compare our theoretical calculations with the recent experimental data by Y. Maeno's group [S. 
Yonezawa et al., Phys. Rev. Lett. 100, 117002 (2008)] and show that there is a good qualitative 
and quantitative agreement between the theory and experimental data. 
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Since a discovery of superconductivity in organic 
(TMTSF) 2 X conductors (X=PF 6 and C10 4 ) [1], their 
physical properties have been intensively studied both 
experimentally and theoretically [2,3]. From the begin- 
ning, it was clear that their superconducting properties 
were unconventional. Indeed, it was found [4] that super- 
conductivity was destroyed by non-magnetic impurities, 
which was recently unequivocally confirmed [5]. In ad- 
dition, it was shown [6] that the conventional for s-wave 
superconductivity Hebel-Slichter peak was absent in the 
NMR experiments. Note that the experimental results [4- 
6] provide strong arguments that superconducting order 
parameter changes its sign on a quasi-one-dimensional 
(Q1D) Fermi surface (FS) of (TMTSF) 2 X superconduc- 
tors. On the other hand, they do not contain information 
about a spin-part of a superconducting order parameter 
and, thus, do not distinguish between singlet and triplet 
pairings. At the moment, the problem about a spin part 
of a superconducting order parameter in (TMTSF) 2 X 
conductors is still controversial. Indeed, early measure- 
ments of the Knight shift in (TMTSF) 2 PF 6 conductor 
[7] in a magnetic field H = 1.43 T showed that spin sus- 
ceptibility was unchanged through the superconducting 
transition. These data were interpreted as an evidence 
for triplet superconductivity. On the other hand, more 
recent Knight shift data in (TMTSF) 2 C10 4 conductor 
[8], obtained in a magnetic field H = 0.957 T, were in- 
terpreted in favor of a singlet superconducting pairing. 

Another unconventional feature of superconductivity 
in (TMTSF) 2 X conductors is very large upper critical 
fields for a magnetic field parallel to their conducting 
planes and perpendicular to their conducting chains, 
H || b' [9-14]. These fields exceed both the quasi-classical 
orbital upper critical field [15,16] and so-called Clogston 
paramagnetic limit [17]. Note that in (TMTSF) 2 PF 6 
superconductor H^ 2 is very large [10,11] due to a for- 
mation of domain walls in the vicinity of antifcrromag- 
netic phase. In contrast, in (TMTSF) 2 C104 conduc- 
tor large upper critical field H^ 2 [12-14] is prescribed 
to 3-D — > 2D dimensional crossover in a magnetic field, 



predicted in Ref.[18] and elaborated in Refs. [19-23]. 
In addition, recent measurements of the upper critical 
fields in (TMTSF) 2 C10 4 superconductor [13,14] have re- 
vealed another unusual property - a novel superconduct- 
ing phase, which appears when a magnetic field is ap- 
plied along the conducting chains, H || a. It is shown 
[13,14] that the above mentioned phase is very sensitive 
to impurities and inclinations of a magnetic field from a 
axis. The authors of the experiments [13,14] have related 
this new phase with a possible formation of the Larkin- 
Ovchinnikov-Fulde-Ferrell (LOFF) state [24,25]. 

The goal of our Letter is to show theoretically that the 
LOFF phase has to appear in a magnetic field, parallel 
to the conducting chains of (TMTSF) 2 C104 supercon- 
ductor, if we use experimentally measured values of its 
Q1D band parameters. The distinctive feature of our 
work is that we take into account both the paramagnetic 
spin-splitting and orbital destructive effects against su- 
perconductivity. This problem is a challenging one, since 
the orbital effects for a magnetic field along the con- 
ducting chains correspond to very particular open elec- 
tron trajectories. These effects cannot be described by 
any of the existing theories, including Refs. [15-23,26]. 
To describe the orbital and paramagnetic effects, below 
we derive an integral equation for a superconducting or- 
der parameter, which, to the best of our knowledge, has 
not been considered before. By means of the measured 
Ginzburg-Landau slops of an anisotropic upper critical 
field in (TMTSF) 2 C10 4 [13,14] and the measured ratio 
t a /tb [27,28], we determine its Q1D band parameters. We 
use these band parameters for a numerical solution of the 
above mentioned integral equation and conclude that the 
LOFF phase has to exist in (TMTSF) 2 C10 4 conductor. 
This conclusion is based both on theoretical arguments 
and on good qualitative and quantitative agreement be- 
tween the theory and experimental measurements [13,14] . 
It is important that the results of our Letter support sin- 
glet d-wave scenario of superconductivity in (TMTSF) 2 X 
materials. 

Below, we consider a Q1D conductor with the following 
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electron spectrum, 

e(p) = -2t a cos(p x a/2) - 2t b cos(p y b*) - 2t c cos(p z c*), 

(1) 

in a magnetic field, parallel to its conducting chains, H j 
a, 



H= (17, 0,0), A = (0,0, Hy), 



(2) 



where t a 3> t b 3> t c correspond to electron hoping inte- 
grals along a , b, and c axes, respectively. 

We represent electron wave functions with definite en- 
ergy and momentum p x in the following way: 

yf(Px]x,y) = exp(±ip x x)exp[±ip±(p x )y] ipf{p x ,y), 

(3) 

where +(-) stands for left (right) sheet of a Q1D FS and 
the functions Py{p x ) are defined by the equations: 



v f (PxTPf) T 2t b cos[py(p x )b*} = 0, 



(4) 



where vf and pp are the Fermi velocity and Fermi mo- 
mentum, respectively. In this case, we can rewrite Eq. 
(1) as: 



Note that finite temperature Green functions for the 
wave functions (8), (3) can be determined by the stan- 
dard equation: 
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(10) 

where uj n is the Matsubara frequency [29] . 

Below, we consider a singlet d-wave scenario of super- 
conductivity in (TMTSF)2C104 conductor, which is in 
agreement with the experiments [4-6,8]. For this pur- 
pose we introduce the following d-wave like LOFF order 
parameter, 

\{x,V,Py) = V2cos(pyb*)exp(iqx)A q (y), (11) 

where the Cooper pairs are characterized by non-zero to- 
tal momentum along a axis, g / 0, and A q (y) is the 
Ginzburg-Landau (GL) order parameter, which depends 
on coordinate y. To derive the gap equation for A q (y), we 
use the Gor'kov equations for non-uniform superconduc- 
tivity [30,31], as it is done, for example, in Ref.[32]. As a 
result of lengthy but rather straightforward calculations, 
we obtain: 



de ± (p) = ±2t b b*(p v -p±)sm(p±b*) -2t c cos(p z c*), (5) 

where we linearize the electron spectrum (1) near the FS 
with respect to momentum p y ; energy 5e = e — ep is 
counted from the Fermi energy ep- 

In a magnetic field, we use the Peierls substitution 
method [18] for Eq.(5), 
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where factor /? takes into account possible decrease of 
the spin-splitting paramagnetic effects due to small de- 
viations from a weak coupling scenario; g is an effective 
electron coupling constant, O is a cutoff energy. Note 
that the last term in Eq.(12) describes the orbital ef- 
fects against superconducting at low enough magnetic 
[TiVyfp^)— -2t c cos( Pz c* - —y) - p B oH\ ipt(p x , y,p z ) nclds ' At hi S h ma S n ctic fields, quantum effects due to 

(•/// Vf [ ...... 



(6) 

where e is the electron charge and c is the velocity of 
light. As a result, we obtain the following Schrodinger 
equation for the electron wave functions: 



= 5e ^t{p x ,y,p z ), 



the Bragg reflections of electrons from boundaries of the 
(7)Brillouin zone, as shown in Refs. [18,32], can improve su- 
perconductivity. This improvement may cause to the ap- 
with ct/2 being a projection of an electron spin on a axis; pearance of the Reentrant Superconducting (RS) phase 
[i B is the Bohr magneton, v y (p y ) = 2t b b* sm(p y b*). [18-23,32]. Our analysis of Eq.(12) shows that the RS 

It is important that Eq.(7) can be exactly solved: phase may appear only at magnetic fields of the order 

.)f 50 T, where superconductivity is already destroyed by 
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he paramagnetic spin-splitting effects. Therefore, for a 
.magnetic field parallel to the conducting axis, H || a, un- 
like the situation, where it is perpendicular to it, H j b' 
j|18-23], we can ignore the above mentioned RS effects. 
From mathematical point of view, this means that we can 



where the wave functions (8) are normalized by condition: 



replace the term sin[ 
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in Eq.(12). If we do this and if we introduce more conve- 
nient variables, y = y, y\ = y + Vy< f^ z, we can rewrite 
Eq.(12) as: 
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It is possible to show that in the vicinity of transition 
temperature, (T c — T) <C T c , Eq.(13) is equivalent to the 
GL expression for the upper critical field, 
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(14) 



where Q{x) is the Riemann zeta function. 

Let us briefly discuss how to determine band param- 
eters t a ,tb, and t c from analysis of the available exper- 
imental data in (TMTSF) 2 C10 4 conductor [33]. First 
of all, the ratio t a /tb ~ 10 was firmly established in 
Ref. [28] , where the measured Lee-Naughton-Lebed oscil- 
lations were compared with the corresponding theoretical 
results. Second, in Refs. [13,14], the GL slopes for the up- 
per critical fields along b and c axes were carefully mea- 
sured: dH b c2 /dT = 2.75 T/K and dH c c2 jdT = 0.14 T/K. 
If we take into account that for d-wave like order param- 
eter (11) in the GL region [33]: 
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(15) 



(16) 



then we can evaluate all three band parameters: t a = 
1340 K, t b = 134 K, t c = 2.6 K. In Eqs.(15),(16), we 
use the following known values: T c = 1.45 K, a = 7.1 A, 
b* = 7.2 A, and c* = 13.1 A [13,14]. 

In Fig.(l), we compare the result of numerical solu- 
tions of Eq.(13) for (3 = 0.9 with the experimental data 
[13] and predict the appearance of the LOFF phase at 
H > 3 T (see also the figure caption). Note that the 
only fitting parameter in our theory is (3. We interpret 
the above mentioned value (3 = 0.9, which fits better the 
experimental data [13], as an evidence of small devia- 
tions from a weak coupling scenario of superconductivity, 
where /J = 1. 

Here, we summarize the main results of the Letter. We 
have derived an integral equation, which allows to de- 
scribe superconductivity in a magnetic field, parallel to 
conducting axes of a Q1D conductor. By using numeri- 
cal solutions of the above mentioned equation for exper- 
imentally measured band parameters of (TMTSF)2C104 
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FIG. 1: Numerical solution of Eq.(13) (solid line) is com- 
pared to the experimental data [13] (dottes). Note that a 
phase transition from uniform superconductivity (q=0) to the 
Larkin-Ovchinnikov-Felde-Ferrell phase (q ^ 0) is predicted 
by our theory in a magnetic field H ~ 3 T. We pay atten- 
tion on a good overall qualitative and quantitative agreement 
between the theory and experiment. 



compound, we have found that the LOFF phase has to 
appear in this superconductor at high enough magnetic 
fields. Comparison of our theoretical curve with the re- 
cent experimental data [13] demonstrates a good overall 
qualitative and quantitative agreement (see Fig.l). 

For theoretical discussions of a possibility of the ap- 
pearance of the LOFF phase in a Q1D conductor in a 
magnetic field, perpendicular to conducting axes and par- 
allel to conducting planes, see Refs. [18,20-23]. There are 
some experimental data in a favor of the existence of 
the LOFF phase in quasi-two-dimensional (Q2D) super- 
conductors k-(ET) 2 Cu(NCS) 2 [34,35], A-(BETS) 2 GaCl 4 
[36], and A-(BETS) 2 FeCl 4 [37]. Nevertheless, theoreti- 
cal calculations, which take into account both the orbital 
and paramagnetic effects in a Q2D case, have not been 
done yet . The paramagnetic spin-splitting effects in the 
LOFF phase in a pure ID case were theoretically studied 
in Ref. [38] , whereas they were theoretically studied in a 
pure 2D case in a number of papers [39-43] . 
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